Properties of solar polar coronal plumes constrained by the 
Ultraviolet Coronagraph Spectrometer data 
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ABSTRACT 

We investigate the plasma dynamics (outflow speed and turbulence) inside polar plumes. 
We compare line profiles (mainly of VI) observed by the UVCS instrument on SOHO at 
the minimum of solar cycle 22-23 with model calculations. We consider Maxwellian velocity 
distributions with different widths in plume and inter-plume regions. Electron densities are 
assumed to be enhanced in plumes and to approach inter-plume values with increasing height. 
Different combinations of the outflow and turbulence velocity in the plume regions are considered. 
We compute line profiles and total intensities of the H I Lya and the VI doublets. The observed 
profile shapes and intensities are reproduced best by a small solar wind speed at low altitudes 
in plumes that increases with height to reach ambient inter-plume values above roughly 3-4 
Rq combined with a similar variation of the width of the velocity distribution of the scattering 
atoms/ions. We also find that plumes very close to the pole give narrow profiles at heights 
above 2.5 Rq, which are not observed. This suggests a tendency for plumes to be located 
away from the pole. We find that the inclusion of plumes in the model computations provides 
an improved correspondence with the observations and confirms previous results showing that 
published UVCS observations in polar coronal holes can be roughly reproduced without the 
need for large temperature anisotropy. The latitude distributions of plumes and magnetic flux 
distributions are studied by analyzing data from different instruments on SOHO and with SOLIS. 

Subject headings: Sun: corona - Sun: magnetic fields - Sun: solar wind - Line: profiles - Plasmas - 
Scattering - Turbulence 



1. Introduction 

This paper is a forward modeling investigation 
of the plasma dynamics of polar plumes. These so- 
called "polar rays" or polar plumes (e.g., van de 
Hulst 1950; Saito 1958; Harvey 1965; Koutchmy, 
1977; DeForest et al. 1997) trace open magnetic 
field lines to different heights depending on ob- 
servation type. In Extreme UltraViolet (EUV) 
emission line spectroheliograms, plumes appear as 
short spikes near the polar limb (Bohlin et al. 
1975a; Ahmad & Withbroe 1977; Widing & Feld- 
man 1992; Walker et al. 1993) with unchanged 
shape and sudden brightness changes over typical 
lifetimes ranging from several hours to days (see 
DeForest et al. 1997; Del Zanna et al. 1997). 



EUV observations can be used to characterize the 
bases of polar plumes and also their plasma char- 
acteristics (densities, temperatures, flows, abun- 
dances, etc.) and also other processes such as 
heating, wave propagation, etc. (e.g. Bohlin et 
al. 1975b; Walker et al. 1988, 1993). Soft X-rays 
show mainly the hot plasma at the foot-points of 
the plumes (e.g., Ahmad & Webb, 1978) that are 
also seen as weak radio sources (Gopalswamy et 
al., 1992). DeForest et al. (1997, 2001a) have 
shown that most white light polar rays can be as- 
sociated directly with the EUV polar plumes and 
their foot-points on the solar surface. 

Above the solar limb, polar plumes appear 
brighter than the surrounding media, suggesting 
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that they are denser than the background corona 
or inter-plume regions (e.g. van de Hulst 1950; 
Saito 1958 & 1965; etc). The coronal-background 
electron density is about ~ 10 8 cm~ 3 and falls 
sharply with height (see different density models 
in Raouafi & Solanki 2006). In comparison, the 
density of polar plumes ranges from 10 8 — 10 10 
cur 3 (Del Zanna et al. 1997; Young et al. 1999). 
Young et al. (1999) reported electron densities 
corresponding to the strong brightcnings at plume 
foot-points of the order of 2.5 — 5.6 10 9 cm~ 3 
and found no decrease with height up to 70 Mm. 
They also inferred temperatures of the order of 2 
MK. In a recent publication, Wilhelm (2006) us- 
ing off-limb spectroscopic measurements from the 
spectrometer SUMER (Wilhelm et al. 1995) on 
board SOHO (Domingo et al. 1995) found that 
plumes are approximately 5-6 times denser than 
inter-plume regions up to k, 1.2—1.3 Rq. The 
plumc-intcrplume density ratio changes little in 
this height range. Note that these are off-limb 
measurements and they do not take into account 
the plume area projected on the solar disk. 

The contribution of polar plumes to the fast so- 
lar wind arising from the polar coronal holes has 
been a subject of debate and controversy. In gen- 
eral, plumes are thought to harbor smaller out- 
flow speeds than inter-plume regions (e.g., Hab- 
bal 1992; Wang 1994; Hassler et al. 1999; etc). 
Thus, polar plumes observed above the solar limb 
seem to be nearly in hydrostatic equilibrium and 
show small or no line shifts, suggesting that the 
bulk of the acceleration of the solar wind takes 
place at larger heights (see Del Zanna et al. 1997; 
Kohl et al. 1997; Wilhelm et al. 2000; Tcriaca 
et al. 2003). Dopplcr dimming techniques ap- 
plied to SOHO/UVCS (Kohl et al. 1995) obser- 
vations of the O VI lines have shown outflows in 
inter-plume dark lanes ranging from 100 to 150 
km s _1 at 1.7 Rq from Sun center and smaller 
speeds in bright plumes (0-65 km s _1 ; see Gior- 
dano et al. 2000). However, Gabriel et al. (2003) 
used the Dopplcr dimming technique applied to 
SOHO/SUMER data and found the intensity ra- 
tio of the O VI doublet to be smaller in plumes 
than in inter-plumes. They concluded that out- 
flow speeds of the solar wind are greater in plumes 
than in inter-plumes. Unfortunately, they did not 
analyze line profiles, in particular the line widths, 
on which the Doppler dimming effect depends. 



EUV coronal emissions in spectral lines that are 
formed within narrow temperature ranges are ex- 
cellent temperature diagnostics for different coro- 
nal structures. SOHO/EIT (Delaboudiniere et al. 
1995) images (Fe XII 195 A that forms at - 1.6 
MK and Fc IX/X 171 A at - 1.0 MK) reveal that 
plumes are cooler than the background corona 
(DeForest et al. 1997; Young et al. 1999). Spec- 
tral analysis have also shown that UV lines are 
narrower in plumes than in inter-plume regions 
(e.g. O VI 1032 A line widths are 10-15% lower; 
Hassler et al. 1997; etc). This suggests that the 
plume material is relatively cooler than or has less 
turbulence than dark lane plasma (inter-plume 
regions). Wilhelm et al. (1998) used intensity 
ratios of temperature sensitive lines for accurate 
measurements of temperatures in the low corona. 
They found evidence for an increasing tempera- 
ture with height in the background corona, but no 
rise of temperature in plumes (see also Young et 
al. 1999). 

The geometry of plumes is an important key in 
understanding their contribution to the fast solar 
wind. The high-speed solar wind at even moder- 
ate solar latitudes is thought to originate within 
coronal holes at high latitudes and to be affected 
by superradial expansion of the coronal hole mag- 
netic field (Zirker 1977). From white light obser- 
vations plumes are thought to also expand super- 
radially with altitude (Ahmad & Withbroe 1977; 
Munro & Jackson 1977; Fisher & Guhathakurta 
1995). DeForest et al. (1997) found that plumes 
expand highly superradially (with a half-cone an- 
gle of 45°) in their lowest 20-30 Mm, and more 
slowly above. DeForest et al. (2001b) also found 
that below ~ 5 Rq plumes expand superradially 
at the same rate as the background coronal hole. 
However, Woo & Habbal (1999) and Woo et al. 
(1999) have claimed to observe radial, rather than 
superradial, expansion of plumes in the coronal 
holes. Speculation on interaction between plume 
and inter-plume material has been reported, in 
particular at high altitudes in the corona. 

Saito and Tanaka (1957a,b) were the first to as- 
sociate plumes with magnetic flux concentrations 
(faculae) in the polar regions. Polar plumes arise 
from footpoints 2-4 Mm (rs 4") wide and expand 
rapidly with height (see DeForest et al. 1997). 
A simulation has shown that a typical plume 
width at the base of the corona is 20 to 28 Mm 
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(Klaus Wilhclm, private communication). Several 
authors (Harvey 1965; Newkirk & Harvey 1968; 
Lindblom 1990; Allen 1994; etc.) have suggested 
that polar plumes arise mainly from unipolar mag- 
netic structures coinciding with the supergranu- 
lar network boundaries. However, DeForest et al. 
(1997) have shown that not all of the flux con- 
centrations give rise to polar plumes. Other au- 
thors suggest that unipolar magnetic flux concen- 
trations reconnect with ephemeral emergent flux 
having opposite polarity (Saito 1965; Wang 1994, 
1998; Wang & Sheeley 1995; etc). This is sup- 
ported by the high temperature observed at the 
plume footpoints, producing high pressure plasma 
from the dipolc that can escape along open field 
lines. DeForest et al. (2001b) found that polar 
plumes are of episodic nature: "polar plumes are 
both transient and persistent: they are recurring 
structures that brighten for only 1 day but reap- 
pear at approximately the same location for up 
to 2-3 weeks" . They speculate that this behavior 
is the result of multiple magnetic reconnections 
of the unipolar flux concentrations at the net- 
work boundaries encountering multiple ephemeral 
dipoles driven by supergranular motion. 

We consider model plumes having different out- 
flow speeds and turbulence and calculate the line 
profile shape of the plume plasma. We use the 
observed profile shapes of the O VI lines as pri- 
mary criteria to select the appropriate height vari- 
ation of those quantities. We then refine them 
to reproduce the other measured quantities, i.e. 
line widths, total intensities and intensity ratios. 
In the following Sections, we study the effect of 
polar plume locations along the line of sight on 
the line profiles and intensities of the O VI dou- 
blet and Lya. The comparison with the observed 
quantities deduced from SOHO/UVCS spectra al- 
lows us to constrain the dynamic properties of the 
plasma inside plumes. In section 2, we present a 
summary of plume spectroscopic observations by 
SOHO/UVCS and highlight the spectral proper- 
ties of O VI 1032 A. A polar plume model is de- 
scribed in section 3 and the results from this model 
are presented and discussed in sections 4-7. Con- 
clusions and prospects are to be found in section 
8. A preliminary version of some of the results has 
been reported by Raouafi et al. (2006). 



2. Summary of UVCS plume observations 

The UltraViolet Coronagraph Spectrometer 
(UVCS; Kohl et al. 1995) has provided impor- 
tant spectral information on coronal lines in polar 
holes. The main spectral lines recorded by UVCS 
are Lya, the O VI doublet at 103 nm and the 
Mg X doublet centered at around 61 nm (in the 
order of decreasing emission). 

Spectral analysis of UVCS data has shown that 
H I Lya has, to a good approximation, a Gaus- 
sian shape at all heights accessible to this instru- 
ment. Kohl et al. (1999) pointed out small devia- 
tions from Gaussians at the profile peaks, allowing 
for two line components resulting from the contri- 
butions of different plasma media with different 
properties along the LOS. Giordano et al. (1997) 
have also shown that Lya profiles are narrower 
when ray paths intersect high density structures 
(plumes) than otherwise (inter-plumes). 

The change with height of the O VI line profiles 
is more complex than for Lya. While at low alti- 
tudes (say below ~ 1.35 Rq) the observed profiles 
are relatively Gaussian shaped, with small devi- 
ations in the wings and peaks, more pronounced 
deviations of the observed profiles from a single 
Gaussian are clearly seen in the observations at 
higher altitudes up to about 2.0 Rq. At this 
height, the O VI profiles are well represented by 
two Gaussian components: a narrow and a broad 
one. The latter appears to be associated with re- 
gions of low density (inter-plume regions) while 
the narrow component shows a clear association 
with the denser regions along the LOS (higher 
emission regions: e.g., polar plumes). The spec- 
tral behavior of the two components as a function 
of height above the limb is quite different. The 
width of the broad component increases dramat- 
ically with height, while that of the narrow one 
remains almost constant (see Kohl et al. 1997; 
1999). Mg X lines have a narrow component at 
r = 1.34 Rq, accounting for only a small fraction 
of the observed spectral radiance. For higher al- 
titudes, the observations are not good enough to 
distinguish whether or not a narrow coronal com- 
ponent exists. Kohl et al. (1999) have shown also 
that the Vi / e values for the narrow components of 
O VI and Mg X are smaller than those for H I 
1216 A at the same heights, and the narrow Mg X 
value is smaller than that of O VI. 
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At the solar activity minimum the polar coro- 
nal holes are well developed and bordering stream- 
ers appear not to reach the UVCS observation al- 
titudes, making it very likely that narrow com- 
ponents are contributions of fine structures inside 
the polar coronal hole. Polar plumes are the most 
prominent features at the heights at which UVCS 
observes, which makes them the main source of the 
narrow component observed in the spectral lines. 

Figs. [1] and display spectral profiles of O VI 
1032 A recorded directly above the pole on 1997 
September 19-23 and 1996 June 03, respectively, 
using the ps 40' long UVCS slit. EIT images 
and Kitt Peak Vacuum Telescope magnetograms 
clearly show the presence of polar plumes against 
the solar disk and magnetic concentrations with 
footpoints at 10° - 15° (Fig.HJ and » 20° (Fig.© 
away from the north pole. The profiles in the first 
figure are formed of two distinct components (nar- 
row and broad). Two Gaussians are needed to 
give satisfactory fits to the profiles at r — 1.75 Rq 
and partly also around r = 2.0 Rq, although with 
a small contribution of the narrow component. 
However, the profiles in the second figure are very 
close to a Gaussian in shape, although they show 
a small degree of asymmetry. Thus, the spectral 
profiles depend relatively strongly on the distance 
of the plume footpoints to the pole and also on 
to the brightness of the plume. In other words, 
they depend on where the plume crosses the LOS 
at a given altitude. This difference in the profile 
shapes is a valuable tool to study plasma dynamic 
in these coronal fine structures. 

The characteristics of the plumes' spectral con- 
tribution (narrow component) to O VI lines as ob- 
served by UVCS are summarized as follows: 

• Dominates the profiles below ~ 2.0 Rq and 
decreases above; 

• Not seen in the profiles beyond ~ 2.5 Rq. 
Only the broad component remains above 
this altitude. Note that the noise level in the 
profiles increases with height, which might 
screen small contributions from fine coronal 
structures along the LOS ; 

• Mostly located at the center of the broad 
component which is close to the rest wave- 
length of the spectral lines. Thus, no sig- 
nificant Doppler shifts are observed for the 



narrow component; 

• The width of narrow component hardly 
changes as function of height in the range 
where this latter is observed (see Kohl et al. 
1997). 

3. A Model of polar plumes 

In order to study the influence of polar plumes 
on the LOS-integrated profiles of O VI and H I, 
we consider a number of plumes aligned along the 
LOS, with foot points placed at different angular 
locations in the polar coronal hole (see Fig. [3]). 
The importance of the angular locations of the 
plume footpoints is enhanced because of the super- 
radial expansion of the magnetic field. The foot- 
point location determines the LOS outflow speed 
of the plume plasma at the location where the 
plume crosses the LOS. Based on the work of De- 
Forest et al. (2001b) and others we assume that 
plumes expand superradially at the same rate as 
the surrounding background corona (inter-plumes) 
according to the global magnetic model by Ba- 
naszkiewicz et al. (1998). The electron den- 
sity is considered to be greater in plumes than in 
inter-plumes, following Wilhclm et al. (1998) and 
Wilhelm (2006). The density contrast is largest 
just above the solar surface and decreases gen- 
tly with height, such that it disappears above 
« 7 Rq (see set in Fig. |4|) . The inter-plume elec- 
tron density stratification is obtained from SOHO 
data (SUMER, UVCS and LASCO (Brueckner 
et al. 1995); see Doyle et al. 1999a, b: solid 
line in Fig. |4|). The density profile of plumes 
(dashed line in Fig. [4j is obtained by consider- 
ing a higher density at the footpoints at the solar 
surface (« 9 x 10 8 cm~ 3 ; see Del Zanna et al. 
1997; Young et al. 1999) and a lower tempera- 
ture inside plumes (~ 10 6 K) than in inter-plumes 
(1.5 x 10 8 cm -3 and 1.2 x 10 6 K, respectively, for 
the model by Doyle et al. 1999a, b). The variation 
of both densities as a function of height is given 
by Eq. (1) in Doyle et al. (1999a,b). 

We assume that the considered coronal species 
have simple Maxwellian velocity distributions in- 
side both plume and inter-plume regions, but with 
different velocity turbulence values, a a (r) , that are 
functions only of height in both structures. The 
plume outflow speed, V pp , is assumed to be pro- 
portional to the solar wind speed in inter-plume 
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regions, Vi P , and the proportionality coefficient is 
assumed to depend only on the radial distance to 
Sun center: V pp (r) = g(r) x Vi P (r), where Vi P (r) is 
obtained through the mass-flux conservation equa- 
tion. Because of the lack of conclusive observa- 
tional results concerning g(r) and a StPP (r) we con- 
sider height profiles for each. 

Two cases are considered for the plume velocity 
turbulence: a s ^ pp (r) = constant (« 50 km s _1 for 
the O VI lines following observational values from 
Kohl et al. 1997, 1999) and a height-dependent 
turbulent velocity a s>pp (r) — f(r). The value of 
50 km s _1 is approximately the observed width of 
the narrow component as measured at 1.34 Rq 
(see Fig. 3 of Kohl et al., 1999). In the latter case 
a low value of the turbulence is prescribed at low 
altitudes, which increases to reach the inter-plume 
value above r w 6.0 Rq. Three cases are consid- 
ered for g(r): constant, linear and non-linear. 

4. Test cases of the turbulence and outflow 
evolution inside polar plumes 

The data presented in Figs. Q] and [2] correspond 
to plumes based between w 10° — 20° from the 
pole. In order to compare our simulations with 
these observed profiles, we limit ourselves in this 
section to the very left plume in Fig. [3] that has 
its foot-point centered at 15° from the pole. This 
plume lies between the examples shown in Figs. [1] 
and [2] 

In the present section, we consider the four test 
cases (sections 4.1-4.4) presented in Table 1. Note 
that we have computed the lines also for the re- 
maining combinations of a s , PP and V pp , but noth- 
ing new was to be learned from them. 

For the inter-plume regions, the plasma prop- 
erties are similar to those found by Raouafi & 
Solanki (2004 & 2006). The theoretical formalism 
used here for the line profile computations is dis- 
cussed by Raouafi & Solanki (2004) and references 
therein. In the present section, we limit ourselves 
to discussing the profiles of the O VI 1032 A line. 
O VI 1037 A is also computed and is found to 
behave similarly to the 1032 A line. Lya profiles 
are also computed, but this line is less sensitive to 
the behavior of a s and v and its shape does not 
change as much as that of O VI 1032 A. Thus, in 
contrast to the O VI lines, the Lya profile is close 
to a Gaussian for most plume parameters and at 



most heights. 

4.1. Case 1: Constant velocity turbulence 
and proportional flow 

In a first step, we consider a plasma outflow 
speed inside plumes that is equal to a constant 
fraction of that in the inter-plume regions. Plume 
speeds that are smaller and greater than inter- 
plume values are considered (V pp (r) /Vi P (r) = 
0.25, 0.5, 1.0 & 1.5). 

Qs, PP is height independent 
and equal to 50 km s _1 . 

Fig.[5]displays the profiles of the LOS-integrated 
O VI 1032 A line at 1.7 R Q and 3.0 R Q for differ- 
ent V pp /Vi P values. For small V pp /Vi P , the plume 
contribution (narrow component) dominates the 
profiles at low altitudes, showing some resem- 
blance with the observed profiles (see solid lines 
in Figs. [S^i-b). However, at high altitudes the 
narrow component is still clearly present with a 
strength well above the typical noise level present 
in the observed profiles. For greater values of 
V pp /Vi P (Figs. Efc-d), the plume contribution is 
Dopplcr shifted away from the line peak already 
at low altitudes (due to the location of the con- 
sidered plume relatively far from the pole). This 
produces strongly asymmetric profiles, unlike the 
observed ones. Also, the plume contribution is 
still significant at high altitudes. The present case 
( a s, PP — constant & V pp /Vi P = constant) does not 
reproduce the overall observed profiles. Conse- 
quently, it does not reflect the plasma conditions 
and mechanisms and taking place in polar plumes. 

Nonetheless, the computations suggest that 
outflow speeds in plumes are likely to be much 
smaller (V pp /Vi P < 0.25) than inter-plume values 
at low altitudes. This avoids large contributions 
from plumes far from the central wavelength of 
the observed line irrespective of their footpoint 
location. Plume outflow speeds close to or greater 
than than the inter-plume values at low altitudes 
will be excluded in the rest of the analysis. 

4.2. Case 2: Constant velocity turbulence 
and linear-height-proportional flow 

In this section, we keep a height-independent 
velocity turbulence inside plumes (a SjPP = 50 km s _1 ). 
However, the plume-interplume outflow speed ra- 
tio is now allowed to vary linearly with height, 
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Table 1 

Test cases of the turbulence and outflow height dependence in polar coronal plumes. C, 

C, a AND b ARE CONSTANTS. 



Case a s , P p(r) V pp (r) 

1. C C'x V lp {r) 

2. C 



(ar + b)x V ip (r) if V pp < V ip ; 



Vi p (r) otherwise 

3- f(r) 



(ar + b)x V ip (r) if V pp < V ip ; 



Vip (r) otherwise 
4. f(r) ' g(r)*V tp (r) 



such that 



V„(r) 



(ar + b)x V lp (r) if V pp {r) < V ip (r), 



V ip {r) 



otherwise, 



(1) 



where a and b are constants chosen such that 
(a r+b) is zero or very close to zero at the solar sur- 
face and increases toward unity at high altitudes, 
roughly above 6 Rq. A possible mechanism for 
such an acceleration of the plume plasma could 
be interaction with the surrounding background 
medium having higher outflow speeds. 

We consider the three height profiles of V pp {r) jVipir) 
that are given in Fig. [6^. LOS-integrated profiles 
computed along rays reaching 2.0 Rq and 3.0 Rq 
from disk center are plotted in Figs.|5J>c. The ob- 
tained profiles are to some extent similar to those 
in Fig. [5] for V pp (r)/Vi P (r) = constant. Although 
the profiles at low altitudes now all have the nar- 
row component nearly centered at the broad com- 
ponent, which is close to the observations, those 
obtained at high altitudes have the same signifi- 
cant narrow contributions shifted to the wings of 
the profile. Such a signal is not observed in the 
corona. 

4.3. Case 3: non-linear-height-dependent 
velocity turbulence and linear-height- 
proportional flow 

The plume outflow speed is that of Case 2, 
while the velocity turbulence of the plume plasma 
is assigned a height-dependent function. A varia- 
tion in the velocity turbulence appears natural if 
cool plume and hot inter- plume material interacts. 
Due to the lack of observational results concerning 
the variation of both parameters, we assume, for 



simplicity, similar profiles for the turbulence inside 
plumes and inter-plumes a Si i P (r). Thus, a StPP (r) 
has the following functional form (see Raouafi & 
Solanki 2004, 2006) 



a StPP (r) — k x arctan (r v + fi) + 5, 



(2) 



where K (= 205/7r), 77 (= 1.5), \i (= 4.5) and 
S (= 130) are adjustable parameters whose val- 
ues have been selected in order to get as close as 
possible to the observations. Since there is ob- 
servational evidence for cooler plume plasma than 
in inter-plumes, we assume that a SiPP (r) is smaller 
than a S) j p (r) at low altitudes and converges gently 
towards the large values of a s ^ p (r > 6 Rq) (see 
Fig. [7^). This guarantees the presence of a nar- 
row component at low heights and progressively 
broader plume profiles as we move further above 
the solar limb. 

Figs.[7h-h displays the O VI 1032 A line profiles 
obtained for lincs-of-sight ranging from 1.5 Rq 
to 3.5 Rq above the pole for the plume outflow 
speed profiles shown in Fig.[6H- At altitudes below 
~ 1.5 Rq the contribution of polar plumes to the 
line profile is narrower than that of inter-plumes 
and is approximately centered at the rest wave- 
length of the line. This results in approximately 
Gaussian-shaped line profiles that are slightly nar- 
rower than those obtained without polar plumes 
(see Fig. Eh). For 1.5 Rq < r < 2.5 R Q a profile 
that can be well represented by two Gaussians and 
that resembles the observed profiles is obtained 
(Figs.[7j:-f). The broad component is due to inter- 
plume material, the narrow component to plumes. 
Due to the low plume outflow speeds at these alti- 
tudes, the Doppler shift of the narrow component 
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is low. 

Above sa 2.5 Rq, the increase in the plume ve- 
locity turbulence results in increasingly broader 
plume contributions, so that the line profile (see 
Figs. [7g-h) takes on the approximate shape of a 
single Gaussian for the dashed velocity profile in 
Fig. [5^ (dotted lines in the same figure) . For the 
other velocity models, the plume component is 
narrower, in particular, for the dotted plume ve- 
locity model, so that one can still distinguish two 
line components at high altitudes. Hence, high 
plume outflow speeds at r > 2.5 Rq that are closer 
to inter-plume values fit the observations better 
than smaller ones. 

The present case is in better agreement with 
the observation than cases 1 and 2, although the 
agreement is not perfect, both at very low and 
very high altitudes. E.g. at low altitudes the 
plume contribution does not dominate the LOS- 
integrated profiles. Probably the plume outflow 
speed at these altitudes is still high enough to dim 
the radiative plume component (shifted from res- 
onance due to the small velocity turbulence). Im- 
provement can be obtained by considering a non- 
linear height-dependent V pp (r) /Vi P (r) allowing fol- 
lower and higher plume speeds at low and high 
altitudes, respectively. This is the aim of case 4. 

4.4. Case 4: non-linear-height-dependent 
velocity turbulence and non-linear- 
height-proportional flow: case best 
reproducing observations 

The velocity turbulence in plumes a StPP (r) is 
again described by Eq. @ and is illustrated by the 
triple-dot-dashed curve in Fig. [T^i. The plume out- 
flow speed is now considered to be related in a 
non-linear way to the inter-plume speed: 

U)=ffWxU). ( 3 ) 

where g(r) is illustrated by the solid line in Fig.JHK 
The plume-inter-plume solar wind speed ratio is 
less than 10% below « 2.0 Rq from Sun center, 
and increases abruptly to more than 75% by about 
3.0 Rq, then converges gently toward unity. This 
profile allows more contribution from the plume 
radiative component at r < 2.0 Rq and favors 
broader plume profiles at r > 3.0 Rq. 

The LOS-integrated O VI 1032 A line pro- 
files are displayed in Figs. [8p-h, together with 



the plume (dotted profile) and inter-plume (long- 
dashed line) contributions. The plume contribu- 
tion is now, as expected, more significant at low 
altitudes, in particular below sa 1.7 Rq and dom- 
inates the profile below « 1.4 Rq (see Fig. |Sp). 

The change in the plume contribution is due to 
combined effects of the electron density decrease 
but mostly to the variation of the outflow speed 
and turbulence of the plasma inside plumes. The 
calculated profiles are close to a Gaussian shape at 
r < 1.3 Rq and display increasingly a bi-Gaussian 
shape above this altitude up to r > 2.0 Rq 
(Figs. Et-d). Above 1.3 Rq, the obtained pro- 
files cannot be represented by one Gaussian and 
the two components (narrow and broad) can be 
distinguished even without fitting. The width of 
the plume component increases relatively rapidly 
above « 2.0 Rq and so does its Doppler shift. This 
sudden change in the profile of the plume contri- 
bution allows the two (broad and narrow) compo- 
nents to merge completely above 2.3 Rq. The 
resemblance of the computed profiles to the ob- 
served ones is now uniformly good (compare pro- 
files in Figs. [I]& [8|). Note that greater gas densities 
in plumes above 7.0 Rq than the one considered 
here allow for a greater contribution from plumes 
that results in broader LOS-integrated profiles at 
high altitudes, that are even closer to the observed 
ones. Any remaining small differences between 
computed and observed profiles could also reflect 
the footpoint location of the modeled plume, or 
the possible presence of multiple plumes along the 
LOS in the observations. 

Case 4 best reproduces the observed profile 
shapes. We expect this model to be the closest 
to the plasma conditions in the polar plumes. In 
other words, the plume material remains cooler 
and much slower than the inter-plume plasma up 
to w 2.0 Rq from Sun center, where a rapid in- 
crease in both outflow speed and effective tem- 
perature takes place, leading to plume dynamic 
properties converging towards those of inter-plume 
region within the next solar radius. The rapid 
change in plume plasma properties is likely to 
be due to interaction between the cool and slow 
plume and hot and fast inter-plume materials, al- 
though some authors speculated that this interac- 
tion takes place much higher in the solar corona 
(see introduction). In the next sections, we will 
consider this case only for an extensive study of 
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the profiles and total intensities of the different 
spectral lines considered here. The contribution 
of polar plumes located very close to or on the 
pole will also be discussed. 

5. O VI doublet 

Fig.[9]displays the height dependence above the 
pole of the O VI 1032 A line width (Fig. EK) and 
total intensity (Fig. [9h>) as well as of the inten- 
sity ratio of the O VI doublet (Figs. [3J:-d) from 
slightly above the solar limb to 3.5 Rq from Sun 
center. The profiles underlying this figure were 
computed according to case 4 (section 4.4). They 
were fit with two Gaussians below 2.0 Rq and with 
one Gaussian above. The widths of the broad and 
narrow components and of the single Gaussian fit 
are represented by different symbols. Also plotted 
are values obtained from UVCS spectra together 
with their corresponding error bars according to 
Kohl et al. (1997) and Cranmer et al. (1999). 
The narrow component widths hardly change with 
height up to 2.0 Rq and fit the measured values 
relatively well. The widths of the broad compo- 
nent, however, change rapidly with height, leading 
to a non-Gaussian shape of the composite profiles. 
Higher up the obtained line widths are compara- 
ble to those obtained by Raouafi & Solanki (2006; 
Fig. 8) although they lie somewhat close to the 
measurements, in particular, around 2.0 Rq. 

Total intensities are also comparable to those 
computed by Raouafi & Solanki (2006) although 
the intensities below 2.0 Rq fit the observed val- 
ues for the O VI 1032 A line better. The change 
in total intensities can be more clearly seen in the 
intensity ratio of the O VI doublet. Above 1.5 Rq, 
the ratios obtained from the composite profiles 
(full circles) are improved compared to the case 
without plumes (open circles), in particular, at low 
altitudes. All the computed ratios now lie inside 
the error bars of UVCS measurements. Note that 
only simple Maxwellian velocity distribution are 
used and no anisotropy in the kinetic temperature 
of the O VI ions is considered. Below 1.5 Rq, 
the intensity ratios with plume contribution taken 
into account are smaller than those without plume 
contribution (see bottom panel of Fig. [9|). In our 
computations, this is due to the large electron den- 
sity numbers in plumes compared to inter-plume 
densities and is not a solar wind effect as reported 



by Gabriel et al. (2003 & 2005). In addition, we 
showed previously that plume wind speeds faster 
or just close to inter-plume speeds give line profiles 
different from those observed in the polar coronal 
holes. Here we demonstrate that the fast solar 
wind originates mainly from inter-plume regions 
and that the plume plasma remains slower and 
cooler up to approximately 2.0 — 3.0 Rq from Sun 
center (see Fig. 9). 

6. HI Lya line 

Fig. [TUk displays the velocity turbulence used 
for the computation of Lya (solid line for inter- 
plume and dotted line for plume). We assume the 
same v pp ^ ya (r) /Vi Pt L ya (r) variation of the plume 
velocity turbulence as in the case of the O VI ions: 
small values at low altitudes and a rapid increase 
toward inter-plume values above w 2.0 Rq. The 
plumc-inter-plume ratio of the outflow speed of 
the hydrogen atoms is assumed to be the same 
as for the O VI ions in the previous section (see 
Fig. [8^) . Note that the hydrogen outflow speed in 
inter-plume regions is lower than that of O VI ions 
(see Raouafi & Solanki 2006). 

The computed line profiles of Lya are shown 
in Fig. ITTJb-i (solid-black lines) together with the 
plume (dot-dashed lines) and inter-plume (dot- 
ted lines) contributions for different rays crossing 
the polar solar axis at different heights. The red 
curves are single Gaussian fits to the computed 
profiles. The high quality of the fits at all heights 
demonstrates how close the computed profiles are 
to a Gaussian shape. However, small deviations 
from the single Gaussian shape are noticeable, in 
particular, at low altitudes in the wings and peaks 
of the profiles. Kohl et al. (1999) observed simi- 
lar deviations and interpreted them as due to the 
contribution of polar plumes. Note that slightly 
changing (ps 10% — 20%) the model parameters 
(i.e. the outflow speed and turbulence of the 
plasma inside plumes) does not change much the 
obtained profiles, which remain Gaussian shaped 
to a good approximation. 

The widths of the computed line profiles are 
displayed in Fig. [TUk as a function of the pro- 
jected heliocentric distance of the corresponding 
rays (full circles). The widths obtained without 
plumes (open circles) are also displayed for com- 
parison, together with the best fit to UVCS ob- 
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servations (solid curve) and the error bars of the 
observed values (vertical lines; see Cranmer et al. 
1999). Line widths obtained in both cases (with 
and without polar plumes) are within the errors 
bars of the measured values. However, the widths 
obtained with plumes taken into account fit the 
observed values better. By optimizing the num- 
ber and position of plumes and their properties, 
an even better match to the observations may be 
obtained. 

Fig. ITOl shows the total intensity of Lya. The 
best fit to the observed intensities along with the 
error bars are also displayed for comparison. At 
low altitudes, intensities obtained with plumes 
taken into account are very similar to those with- 
out plumes and are slightly greater than the mea- 
sured ones. However, at higher altitudes inten- 
sities obtained with plumes provide an improved 
fit to the observations. Generally, the presence of 
plumes along the LOS slightly increases the to- 
tal intensity due to the density enhancement in 
plumes. 

7. On the foot-point locations of polar 
plumes 

Fig. [TTa-d shows the LOS-integrated profiles of 
the O VI 1032 A line with contributions from each 
of the four right-hand plumes in Fig. [3] taken into 
account separately. The contributions from the 
two polar plumes based at colatitude —5° and 0° 
(the narrow component of the line profile) domi- 
nate the profiles up to sw 2.5 Rq and remain clearly 
visible up to w 3.5 Rq (dotted and solid lines) 
which is not the case for the observed profiles. The 
reason for the continuing visibility of the truly po- 
lar plumes up to large altitudes is that the rapid 
acceleration of the solar wind in the plume hardly 
affects the plume profile, since it is almost per- 
pendicular to the LOS. The contribution from the 
plume at 10° is important at low altitudes but rel- 
atively small at high altitudes (dashed lines). This 
suggests that polar plumes preferentially originate 
more than w 10° away from the pole. The contri- 
bution from the plume with footpoint at 18° (dot- 
dashed lines) is very similar to that shown in Fig. [2] 
due to plumes at « 20° from the solar pole. 

Various authors have studied locations of po- 
lar plume foot-points within polar regions. They 
traced plumes arising near the edges of coro- 



nal holes, relatively far away from the solar pole 
(eclipse observations; see Saito 1965). In contrast, 
we did not find any report on plumes very close to 
the polar axis. This case may be very rare for un- 
known physical reasons (for example, bipolar flux 
may not emerge near the pole). 

Polar coronal plumes can be easily identified 
off-limb in white light (LASCO, eclipses) and EUV 
(EIT, TRACE) images. On the solar disk, for in- 
stance by using EIT images, it is also possible 
to trace them back to the solar surface, in par- 
ticular, using images recorded in the wavelength 
bands containing the hot lines 171 A and 195 A 
that are regularly recorded by EIT. X-ray im- 
ages also allow the determination of the locations 
of plume foot-points, that correspond very often 
to bright points reflecting the high temperatures 
there. Once the plume foot-points are identified, 
it is relatively easy to link them to the underlying 
magnetic structures using magnetograms. Data 
recorded daily over several consecutive days or 
better still weeks help to survey the latitudes of 
polar plume foot-points. For this purpose, peri- 
ods of time when the solar axis is tilted toward 
the observer are more useful since the observation 
of a major area of the polar coronal hole including 
the pole is possible. These times correspond to 
March and September of each year. 

EIT images recorded at solar minimum seem 
to support this hypothesis. An average of all daily 
unsigned SOLIS (Synoptic Optical Long-term In- 
vestigations of the Sun) high sensitivity magne- 
tograms recorded in September 2005 (see bottom 
panel of Fig. [TTjl show magnetic flux concentra- 
tions preferentially between w 70° and « 80° lat- 
itude, although we are not at solar minimum yet. 
Another reason why plumes are more likely to have 
their footpoints some degrees away from the pole 
has got to do with the relative area available for 
plumes between 70 and 80° is 3 times larger than 
at latitudes above 80°. However, preliminary anal- 
ysis of SOLIS magnetograms tends not to support 
the area argument. This subject will be addressed 
in detail in a subsequent publication. 

8. Discussion 

The aim of the present study is to obtain in- 
formation on the height evolution of the plasma 
dynamics inside polar plumes. We compared 
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profile shapes of the O VI doublet (1032 and 
1037 A) and the H I Lya spectral lines observed by 
SOHO/UVCS at different altitudes in polar coro- 
nal holes to synthetic profiles that include the con- 
tribution of polar plumes placed at different lo- 
cations along the LOS. Observational constraints 
on these fine coronal structures (densities, tem- 
peratures, geometry, ...) have been taken into ac- 
count following different reports in the literature. 
Since the outflow speed and kinetic temperature 
evolution with height are not conclusively known, 
we consider different height profiles for the out- 
flow speed and plasma turbulence inside plumes 
and study how these choices influence the line pro- 
files of the O VI doublet, which are more sensitive 
to the detailed dynamic structure along the LOS 
than H I Lya. 

By comparing the synthetic profiles to the ob- 
served ones in the polar coronal holes, we de- 
termined plausible height profiles for the outflow 
speed and plasma turbulence in polar plumes with 
respect to the surrounding background corona. 
We find that models with constant (height- 
independent) velocity turbulence in the plumes 
show narrow components with a significant contri- 
bution to the LOS-integrated profile at all heights, 
in particular above ~ 3.0 Rq where there is no 
observational indication of the presence of such a 
narrow component in the observed profiles. This 
is independent of whether the plume outflow speed 
equals a constant or a height-dependent fraction of 
the inter-plume outflow speed. Synthetic profiles 
with shapes very similar to the observed ones are 
obtained when the plasma turbulence in plumes is 
height-dependent with similar variation to that in 
inter-plume regions and converges towards the lat- 
ter at greater heights. At the same time, the out- 
flow speed is small near the solar surface, but ap- 
proaches the inter-plume value at greater heights. 
The narrow component is found to contribute sig- 
nificantly to the resultant profiles at altitudes be- 
low ^2.0 Rq and is increasingly dimmed above 
that altitude. This model suggests that the plume 
plasma remains much slower and cooler (in the 
sense that a SjPP < a Sjip ) than inter-plume ma- 
terial up to w 2.0 Rq. Above « 2.0 Rq these 
parameters experience a rapid rise towards inter- 
plume values reached by w 3.0 Rq. The rapid 
acceleration and increase in the effective temper- 
ature of the plasma inside plumes is probably the 



result of interaction with the faster and hotter 
inter-plume material. The entrainment or mixing 
of material between plumes and inter-plumes is 
hindered by the magnetic field. The smaller the 
magnetic energy density relative to the thermal or 
the kinetic energy density, the more efficiently en- 
trainment takes place. Now consider the following 
numbers: from r = Rq to r = 2i? Q the density p 
decreases by a factor of approximately 10 3 , while 
the area of a plume increases by roughly VlO 3 . 
Consequently the kinetic energy density, pv 2 /2, 
remains roughly constant over this height range, 
while the magnetic energy density, B 2 /2p drops 
by roughly a factor of 10 3 . This may explain why 
the plum properties remain distinct from those 
of the inter-plume gas up to about r — 2Rq and 
then start to merge. As far as the actual interac- 
tion is concerned, it may be enhanced by, e.g., the 
Kelvin- Hclmholtz instability which can act at the 
boundary between two streams with sufficiently 
different speeds. In addition electron collisions 
may not be sufficient at heights > 2.0 Rq to keep 
the different species in equilibrium. Consequently, 
plume species are heated and accelerated differ- 
ently and rapidly as the density drops sharply 
with height. This may not be the case below this 
altitude where the density is high enough to keep 
all species in equilibrium contrasting inter-plumes 
regions where the (electron) density is smaller al- 
lowing for the plasma equilibrium to brake down 
at lower altitudes and thus the hotter and faster 
plasma compared to plume regions. 

We find that observed profile shapes, line 
widths and total intensities of the O VI doublet 
and of Lya are better reproduced when the influ- 
ence of appropriate plumes is included. In par- 
ticular, the intensity ratios of the O VI doublet 
below 1.5 Rq are smaller when the plume contri- 
bution is taken into account. This is due to the 
density enhancement in the plume region. This 
result is in disagreement with the interpretation 
preferred by Gabriel et al. (2003 & 2005) that 
such a reduction in the ratio is due to the Doppler 
dimming effect, which would imply a higher out- 
flow speed in plumes than in inter-plume regions. 
In our analysis, such high outflow speeds in the 
plumes produce synthetic profiles whose narrow 
component is Doppler shifted towards a profile 
wing, in contrast to the observations. 

It is noteworthy that including polar plumes 
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increases the agreement of the computed profiles 
with the observed ones and with critical param- 
eters like line width or O VI line ratio even if 
no anisotropy in the temperature (or turbulence 
velocity) is included in the computations. Note 
that the properties of the plumes were constrained 
basically using the line profile shape only and 
the quantitative comparison with the line param- 
eters and ratio was only carried out later. This 
strengthens the thesis put forward by Raouafi & 
Solanki (2004 & 2006) that there is no compelling 
need for large anisotropies in temperature in coro- 
nal holes based on published UVCS data. 

Note that we have not tried to optimally repro- 
duce the observed quantities by carrying out an 
exhaustive search in parameter space. Thus the 
obtained results can possibly be improved further 
by applying a x 2 minimization procedure to fit the 
line profiles. 

Plumes with footpoints within 10° — 15° of the 
pole display significant narrow components even 
at altitudes above w 2.5 Rq, contrary to the ob- 
servations. EIT images during the solar minimum 
together with high sensitivity magnetograms from 
SOLIS show evidence supporting the hypothesis 
that plumes originate away from the pole. This 
subject will be studied in detail in a subsequent 
paper. 
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Fig. 3. — Spatial distribution of the polar plumes 
used for the calculations of the coronal line profiles 
in the plane of the LOS. The plume footpoints are 
18 Mm wide and are centered at angular positions 
relative to the pole from left to right —15°, —5°, 
0°, 10° and 18°, respectively. The magnetic field 
is assumed to expand similarly in both plume and 
inter-plume regions. No plumes are placed further 
than 20° since these plumes do not reach large 
apparent heights above the pole due to the rapid 
expansion of the magnetic field. However, their 
presence may affect profiles at low altitude (below 
«2.5i? ). 
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Fig. 4. — (a): Electron density stratifications in- 
side plumes (dashed line) and inter-plumes (solid 
line). The inter-plume density is that obtained 
from SOHO (SUMER, UVCS and LASCO) data 
(see Doyle et al. 1999a, b). The plume density 
stratification is obtained by imposing a greater 
density at the solar surface (sec text for refer- 
ence) and lower temperature in plumes relative 
to the surrounding ambient corona, (b): Plume- 
interplumes density ratio as function of height. 
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Fig. 7. — Case 3: (a): a SiPP (r) and a Siip (r) 
as a function of height used for the computa- 
tion of the profiles in frames b-h. (b-h): O VI 
1032 A LOS-integrated profiles corresponding to 
the three plume velocity models shown in Fig. 6 
are given with the same linestyle as in that figure. 
Plume (smaller, narrower profiles at bottom of 
each frame) and inter-plume (long-dashed curves) 
contributions are also plotted. 
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Fig. 1. — Top: EIT 171 A images of 1997 Sep. 19 and 23, respectively. Polar plumes with strong emissions 
are present more than 10° away from the pole on both days (latitude marked every 5°). Middle & bottom: 
O VI 103.2 nm line profiles (+ signs) recorded by UVCS over the same period of time. The profiles are 
better fit by two Gaussians (solid curves) except for the bottom-right panel where a one-Gaussian fit is 
equally good. The narrow and wide fit components are represented by dotted and dashed lines. Counts in 
the different panels are not normalized to exposure times. 




Fig. 2. — EIT 171 A image (top-left) and Kitt Peak vacuum telescope magnetogram (top-middle) of the 
north solar pole on June 3, 1996. Latitude lines indicate 5° increments. Polar plumes with footpoints about 
20° away from the pole are clearly visible in the EIT image. The footpoint of the plumes around pixels 
45 x 10 in the top-left panel lies in the white box in the magnetogram. The other panels display the line 
profile of O VI 103.2 nm recorded by UVCS on the same date at different heights (+ signs). The profiles are 
slightly asymmetric although they are fit relatively well by a single Gaussian (solid curves). 
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Fig. 5. — Case 1: LOS-integrated model profiles of the O VI 1032 A line at 1.7 Rq (solid lines) and 3.0 Rq 
(dotted lines) in the presence of the very left polar plume in Fig. [3] The plume velocity turbulence is assumed 
to be constant and the outflow speed is equal to a height-independent fraction of the inter-plume outflow 
speed as noted in the different panels. 
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Fig. 6. — Case 2: O VI 1032 A line profiles for a wind speed in plumes equal to a linear-height-dependent 
fraction of the inter-plume outflow speed and a constant velocity turbulence in plumes, (a) The chosen flow 
speed height profiles, (b) resulting line profiles at 2.0 Rq and (c) at 3.0 Rq from Sun center. 
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Fig. 8. — Case 4: The same as Fig. but for plume outflow speed equal to a non-linear height-dependent 
fraction of inter-plume speed (a). The line profiles computed along different rays now show considerable 
similarity to the observed ones (b-h). The contribution of plumes (dotted profiles) is more significant at 
low altitudes and is extensively dimmed relative to the inter-plume contribution (dashed profiled) at higher 
altitudes. Line profiles are fit by two Gaussians up to 2.0 R@ and by a single Gaussian for greater heights. 
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Fig. 9. — Model (circles and triangles) line width (a), total intensity (b) and intensity ratio (c,d) of the O VI 
doublet vs. r. Values obtained from synthetic line profiles following case 4 (profiles plotted in Fig 8) are 
shown as filled and open circles and triangles. The solid and dashed lines are best fits to UVCS observations, 
the solid error bars represent the accuracy of the observations (see Cranmer et al. 1999). The dot-dashed 
vertical bars represent the uncertainty in the observed narrow component widths (see Kohl et al. 1997). 
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Fig. 10. — (a): velocity turbulence of hydrogen atoms in a plume (dotted line) and in inter-plume (solid 
line) as a function of height, (b-j): computed Lya profiles: composite profile (balck-solid lines); the plume 
and the inter-plume contributions (dotted and dot-dashed lines, respectively) and single Gaussian fits to the 
composite profiles (gray-dashed lines). Computed profile widths (k) and total intensities (1) of Lya with 
(full circles) and without (open circles; see Raouafi & Solanki 2006) polar plume crossing the LOS. The 
solid curves in (k) and (1) are the best fits to UVCS measured values, together with the error bars of the 
measurements (see Cranmer et al. 1999). 
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Pixel 

Fig. 11. — (a-d): line profiles of the VI 1032 A line at different altitudes with the individual contribu- 
tions of the four right-hand plumes in Fig. [3] taken into account separately. The intensity unit is photons 
s _1 cm~ 2 sr _1 A -1 . Contributions from plumes at or close to the pole are noticeable beyond 2.5 Rq which 
is not the case for the observed profiles. Note the great similarity of the dot-dashed profiles (with contri- 
bution from the farthest plume to the pole (18°)) with the observed profiles in Fig. [2] Bottom: SOLIS 
magnetograms averaged over September 2005 showing the unsigned flux distribution around the north polar 
coronal hole. The flux is concentrated preferentially between 70°and 80°latitudes. The solid curves mark 
latitudes every 5°. 
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